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Using perturbative QCD, the contribution at the leading twist, leading a 3 level, of charm and up quark loops to the decays 
B —* K*^y and B — * fry is presented. In the case of B — * fry, the relative importance of these contributions depend upon the 
unknown CKM matrix elements Vbu and Vtd- Assuming that the ratio r — Vb c V* d /VbtV t * d is bounded between —2.25 < r < —.5 
as is suggested by the Particle Data Group, the error in extracting |Vtd/Vt 3 | by these decays is estimated. 
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I. INTRODUCTION 

The recent observation Q of the rare decay B — > if*7 
was the first unambiguous experimental verification of 
flavor-changing neutral currents. While not occurring 
at tree-level in the Standard Model, such currents, or 
"penguins" , are well-known || to arise due to loop ef- 
fects. The magnitude of these transitions depend || upon 
fundamental parameters in the Standard Model such as 
the mass of the top quark as well as the CKM matrix ele- 
ments VtbV t *, although significant QCD enhancements ||] 
tend to offset somewhat the dependence on the former. 
A systematic determination of the two-body rare decays 
of the B meson, both for determining these fundamental 
parameters as well as for indications of possible physics 
beyond the Standard Model, forms a major part of the 
goals of the future B factories at SLAC and KEK. De- 
spite these connections, significant theoretical input will 
nevertheless be necessary to interpret the data so as to 
distinguish purely hadronic effects as well as potentially 
competing mechanisms. 

With the latest data on B — > J/ipK* ||, the ques- 
tion of the relative importance |jj of "long-distance" ef- 
fects arising from intermediate charmonium states to the 
decay B — > if* 7 has received new attention 0. The 
potential relevance of annihilation diagrams in the decay 
B — > p7 has also been emphasized by Atwood et al. as 
a possible source of contamination of the extraction using 
these two decay modes j| of the ratio |Vtd/Vt a |. While 
these effects are undoubtedly present, it is crucial for 
their analysis that all mechanisms are understood within 
a single formalism. For example, most of the studies 
cited employ vector-meson dominance to estimate from 
the measured charmonium decay rate a contribution to 
the B — ► K*"f transition amplitude purportedly distinct 
from the b —> sj penguin vertex already present in the 
effective Hamiltonian. However the (significant) running 
of the Wilson coefficient of the electromagnetic penguin 
operator in the effective Hamiltonian is predominantly 
driven by its mixing with other operators via charm- 
quark loops. There is hence a serious potential of double- 
counting when using a mixed meson-quark language to 
describe the transition amplitude. This potential is fur- 



ther emphasized when recalling that the mixing between 
7 and O2 (see Eqs. (0) and ( [Tl] ) below) is in fact iden- 
tically zero at the one-loop level. 

Recently, perturbative QCD (pQCD) methods were 
applied to the penguin decay B — > K*j |^0| . Previ- 
ously pQCD had been found |llj to be quite successful 
Jl2| in describing the hadronic, two-body decay chan- 
nels of the B meson. The large mass scale of the de- 
caying B meson, coupled with the restriction to the two 
body exclusive decay modes involving nearly massless, 
highly Lorentz contracted states, implies that the tran- 
sition amplitude is governed by short-distance processes 
(on hadronic scales) and hence an appropriate environ- 
ment to apply pQCD. Indeed such simple kinematic con- 
siderations raises the question as to the true suitability 
of calling any mechanism contributing to these decays 
"long-distance" and further emphasizes the need for a 
unified approach. 

The dominant graph in the pQCD framework that con- 
tributes to the decay B — > K*j is shown in Fig. (1). 
A branching ratio of roughly 3 x 10~ 5 was obtained in 
Ref. |l(J , in reasonable agreement with the data jj] . One 
process that was however omitted in that analysis was the 
so-called double penguin graphs of Fig. (2). Based on the 
analysis of Ali and Greub || to the decay b — ► s 7 g, it was 
argued that such graphs would be suppressed compared 
to the dominant decay mechanism shown in Fig. (1). The 
fact that the virtuality of the gluon in Fig. (2) is in gen- 
eral significantly off-shell (of O(m^) in practice) and is 
embedded in additional loops, implies that the analytical 
expressions found in Ref. || using on-shell gluon emis- 
sion could not be simply applied. A reasonable founda- 
tion for omitting the graph was therefore welcomed and 
exploited. 

For a detailed understanding of the decay amplitude, 
and especially for the extrapolation to the case B — > fry 
where CKM factors no longer suppress other background 
processes, it is nevertheless important to quantify these 
additional mechanisms. This paper reports the results 
of such an analysis. While as we will see, these re- 
sults do in fact support the conjecture that these mech- 
anisms enhance the total decay rate, nontrivial hadronic 
phases enter into obtaining this result. This possibility 
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has been generally ignored. The enhancement is related 
to the observed || failure of factorization in the decay 
B — > J/i/jK* , although the overall enhancement factor of 
the decay rate itself depends sensitively on details of the 
wavefunction of the K* in a fashion that the dominant 
decay mechanism, Fig. (1), does not. 



II. CALCULATIONS 

Exclusive processes at large momentum transfer are 
addressed |l3| within pQCD starting with a Fock com- 
ponent expansion of the involved hadrons whereby a twist 
expansion suggests that the contribution from the low- 
est order Fock component dominates the physical observ- 
able under consideration. An exclusive process then in- 
volves a perturbatively calculable hard amplitude convo- 
luted with a nonperturbative, soft physics wavefunction, 
tjj m , from each of the hadrons m entering or leaving the 
hard interaction. These wavefunctions, although as yet 
uncalculable from first principles, are universal for each 
meson, i.e. they factorize from the hard amplitude and 
hence are independent of the process involved. In exclu- 
sive processes in pQCD they play the analogous role that 
structure functions do in the case of inclusive scattering 
events. Thus as was employed in Ref. ideally one 
can phenomenologically parametrize these wavefunctions 
using a (few) measured cross-sections/decay rates. 

The factorization scheme advocated by Brodsky and 
Lepage Jl3| is employed, whereby the momenta of the 
quarks are taken as some fraction x of the total momen- 
tum of the parent meson weighted by a soft physics dis- 
tribution amplitude 4>(x). The peaking approximation is 
used for (j>B, the distribution amplitude of the B meson, 
wherein 



Jb_ 

2%/3 



S(x - e B ) 



(1) 



The decay constant of the B is fs (in the convention 
/tt = 93MeV) and x is the light cone momentum fraction 
carried by the light quark. The parameter cb in 4>b(x) 
is related to the difference in the masses of the B meson 
and b- quark, 



m B = m b + Ai 



(2) 



whereby es = Ab/ttib- 

Ignoring the mass of the K* , its distribution amplitude 
can be written as 

<t>K»[y) = V3f K *y(l- y)4>K-{y)- (3) 

Two guesses for <pK*(y) were considered in |Tc| ], 

<pK*(y) = l 



The first is the so-called super-asymptotic |L4| distribu- 
tion amplitude, while the lower form is one suggested by 
Chernyak, Zhitnitsky and Zhitnitsky (CZZ) ]l5| for the 
transverse polarizations of the p meson, tl 

In the present context, the factorization scheme is aug- 
mented by the viability of an €b expansion for the decay 
amplitude. All terms in the hard amplitude of order e B 
are ignored, both because they are expected to be nu- 
merically small and because they are related, through 
the mass of the light quark, to transverse momentum 
effects. Since the latter is ignored in the factorization 
scheme, self-consistency dictates these other terms also 
be ignored. The parameter es has been fitted [|ll|jl7| 
using the decay B — > Dir. With mild assumptions con- 
cerning the decay constants Jb and fjj, a typical value 
found was eb = .095. 



A. Dominant mechanism 

In Ref. pel] , Fig. (1) was shown to dominate all graphs 
involving penguin operators. To leading order in e^, the 
contribution to the decay B + — > K* + j was found to be 



&GVtbV£ 
m B e B 



a s {fi)C 7 (p)I 



x {p-qe*-e+ie^ a(i p»q v e* a e ), (5) 



where e (£) is the polarization of the photon (K* 
eC F G F 



(6) 



and Ci{n) is the renormalization group improved Wil- 
son coefficient of the electromagnetic penguin operator 

° 7 = T^ m ^ F ^\ (1 + 75) b. (7) 

The quantity I involves an integral over the distribution 
amplitude of the K* and is given by 

1= / dy$ K *{y) 7. 777 ■ 8 

Jo y~ 2e B - i0+ 

For the two distribution functions considered, this be- 
comes (using e B = .095 as discussed earlier) 



I = 



.68 + i 2.55 
1.75 + i 1.96 ' 



(9) 



where the upper number is for the asymptotic distribu- 
tion and the lower is for the CZZ one. The imaginary 
part comes from an internal propagator going on-shell, 



<P K ,{y) = by 2 (l-yf 



(4) 



'See however Ref. M for a QCD sum rule result for <fi K * {y) 
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kinematically allowed here by the inputed soft-physics 
parametrization Eq. (|^) that mt < M B . As in other 
cases in pQCD Q , it is a calculable hadronic phase since 
the overall kinematics of the reaction dictate that only 
short distance propagation occurs, as discussed in [ fL7| . 
In more technical language, the pole is not pinched and 
hence not associated with a long distance event |ll| . 

One curious result of Eq. ^ is that I 2 is nearly iden- 
tical despite the significant differences in phase for the 
two distribution amplitudes under consideration. Hence 
the total decay rate was found in [|l0| insensitive to these 
soft physics inputs. This pattern however does not con- 
tinue when additional mechanisms are included in which 
interference effects depend crucially on the details of the 
relative phases. 



B. The Annihilation graphs 

The first such competing mechanism, the annihilation 
diagrams of Fig. (3) was already considered in (To). Ill As 
suggested by the figure, the dominant result (in l/e B ) 
is from the graph involving photon emission off of the 
light-quark of the B meson. To this order in e B , the 
contribution of this process to the decay B + — > K* + j is 



2c, 



M a . nn = 



(c 2 (A) + 1 ^c 1 (A) ; mK . 

m B e B m B 

G F 

e^j=V u bV* s f B fK* 
x (p • qe* • + ie^ a ^q"e* a C P ) 



(10) 



where e u = 2/3 and C2 and C\ are the Wilson coeffi- 
cients of the 4-point operators (greek subscripts are color 
indices) 

01 = iu Q 7 M (l - 75)6,3 S/37 M (1 - 7 5 )w Q , 

02 = ^w Q 7 M (l - 7 5 )& Q ^7^(1 - 7 5 )w^. (11) 

In this case vtik * is kept as it appears as an overall fac- 
tor arising from the W turning into a K* . Due to the 
differences in CKM matrix elements as well as the factor 
m^/Afs, the annihilation amplitude is here essentially 
ignorable. As though discussed by Cheng and also At- 
wood et al. this is no longer true for the (as yet un- 
seen) decay mode B — > pj. For the decay B + — ► p + 7, 
M ann is obtained from Eq. ([l0|) with the obvious modi- 
fication in CKM factors and meson decay constant. For 
the neutral decay B° ~ > p° 7, 



^The importance of these annihil atio n diagrams was first 
realized by Bander et al. in Ref. [|L9| in the context of D 
meson decays using a nonrelativistic quark model approach. 



led (Ci(m) + T^^Mm) j m h 
m B e B m B 

e—V uh V ud!B — 
x (p ■ qe* ■ C + ie^ a ^q"e* a C P ) 



(12) 



C. The Quark-Loop graphs 

We now come to the second competing mechanism, the 
quark-loop graphs of Fig. (2). These have been recently 
considered by Greub et al. in Ref. po|] using a quark 
model approach. These authors were predominantly in- 
terested in studying CP violating effects and thus fo- 
cussed on only the absorptive parts of these graphs. Here 
though, the entire amplitude is of interest. The impor- 
tant application to CP violation will be deferred to a later 
work. 

In the evaluation of the graphs of Fig. (2), one finds 
that while each graph is individually ultraviolet diver- 
gent, their sum is finite. Likewise gauge-invariance (in 
both the strong and electromagnetic interactions) is only 
obtained after the graphs are summed. In the case that 
the gluon is on-shell, Q 2 = 0, entirely analytical results 
|}| are possible. Such a convenient form has not been 
found in the present context which involves a second 
integration over the gluon's virtuality, Q 2 — —ye B M B . 
Intermediate results will be presented, in which the re- 
maining integrals were then evaluated numerically. The 
contribution of the two quark loop graphs of Fig. (2), 
M q q loop, to the decay amplitude is 



M„ ff 



qq loop 



16 ^bV* a f dy^ K ,{y)I{y) 

3m B e B J 



x(p-qe*-C + ^ a0 p^e* a C ) 



(13) 



where I{y) is 
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1 ( -ml 

dx — In 

\2q-Qx 



m 2 + Q 2 x(x - 1) 



+ (1 - x) In 
2q ■ Qx 2 



m 2 + 2q ■ Qx(x — 1) 
m 2 + 2q ■ Qx(x - 1) 



In 



m 2 + Q 2 
(1 - z+)(x - Z-) 



1 



(1 - zJ)(x - Z+) 

2" 

- iirQ(y - y a ) x 




Vl-4a m x+ 
2 2q ■ Q x 



x T 
dX J_ 



(14) 
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in which 

Vo -- 



4m^ 
TFT 



z± = 



_ -2q-Qx±/3 



2Q 2 



X ± = 



M 2 B [y 

1 ± VI - 4a 



(15) 



and /3 = w (2g • Qs) 2 — 4Q 2 m 2 . To obtain this result, it 

is important to strictly keep only terms that are leading 
in e^ 1 . Subleading terms are higher-twist and require 
other higher-twist elements (such as transverse momen- 
tum degrees of freedom) to maintain gauge-invariance. 
Numerical integration yields that 



dy4>K*(y) l(v) 



.51 
.43 



i .14 
i .61 



.71- 
.64 



« .30 
i.30 



(16) 



where the results are presented as in Eq. (g). The two 
columns are the results for the case of the up (m q = 0) 
and charm quark loops (m c = 1.5GeV) respectively. 
There is negligible change in the charm quark results tak- 
ing l.25GeV <m c < 1.75GeV. 



III. RESULTS 



For the decay B — > K* 7, the two relevant ampli- 
tudes are -/V/^peng and Afccioop- As indicated by their 
expressions, in order to reflect the fact that the aver- 
age virtuality of the exchanged gluon (Q 2 ) in each of 
the mechanisms is different, the running coupling and 
Wilson coefficient are evaluated at scales appropriate 
to each amplitude. For the dominant piece, M 7peng , 
this occurs at p 2 s» .5GeV 2 ; — Q 2 at the pole in the 
bottom quark's propagator. In the case of M c gi oop , 
(-Q 2 ) ~ 1^bM% = 1.3GeV 2 . Note that a change in 
scale /i 2 affects the overall magnitude of the branching 
rate predominantly by the dependence on a s in Eq. (^|). 

For Aqcd — -2GeV one obtains for the enhancement 
factor R in the total decay rate the result that 



R 



(A^7p en g -M cc loop) 

M 2 

7pcng 



1.06 
1.31 



(17) 



where the two results are for the asymptotic and CZZ 
distribution amplitudes respectively. Note that in this 
ratio the uncertainty due to our present ignorance of fs , 
as well as the dominant dependence on the parameter cb 
cancel. There is little dependence in R on the exact value 
of Aqcd- However a precise value for R is clearly depen- 
dent upon details of the kaon's soft-physics information 
(unlike the square modulus of Eq. ([)])). The range is due 
primarily to the phase differences in .M 7P eng (Eq. Jl6| ) 
shows that M cE i oop is nearly insensitive to the choice of 
(f>K*{y))- One should note that this range is compara- 
ble to the various estimates obtained using vector-meson 



dominance methods |6J,|7|. Such duality is perhaps best 
understood by further noting that in pQCD [^l| the non- 
factorizing amplitudes are found to produce large contri- 
butions to the transversally polarized final states of the 
decay J/ip K W. 

In the case B + — > p + 7, CKM factors no longer sup- 
press either the annihilation diagrams or the up-quark 
loop contributions. To quantify this dependence, unitar- 
ity of CKM matrix is exploited, whereby 



vm + v bc v; d + v bt v t 



'td 



0. 



(18) 



Defining the ratio Vb c V* d /VbtV t * d = r and assuming SU(3) 
symmetric distribution amplitudes for the p and K* 
mesons one obtains for the ratio of decay rates that 



B+- 



fpVtd 

fx-Vu 



fa 



sy.CZZ 



(r), 



(19) 



where f a sy,czz{f) 7^ 1 represents the error due to com- 
peting mechanisms of extracting this ratio. The sub- 
scripts recall the fact that f(r) is dependent upon the 
distribution amplitude used for the p and K* mesons. 
In Fig. (4), f(r = -2.25) and f(r = -.5) are plotted 
as a function of Aqcd (the range of r as suggested by 
the Particle Data Group is -2.25 < r < -.5, |§|]). We 
observe that while f(r) is fairly independent of the dis- 
tribution amplitudes, there is significant dependence on 
Aqcd- Note that the deviation of f(r) from 1 is greatest 
for smallest Aqc d , where one would expect that the per- 
turbative formalism employed here is more accurate. A 
conservative estimate is therefore that these these decay 
modes can be used to obtain |Vtd/Vt a | || to only within 
a factor of two or so. 

In the case of the neutral B decays, the effective ab- 
sence of the annihilation diagram improves upon this re- 
sult. From Fig. (5) we obtain an estimated uncertainty 
in extracting |V t d/Vt s | on the order of 50%. 



IV. CONCLUSION 

The two body exclusive decays of the B meson nec- 
essarily involves modeling of the initial and final state 
hadrons. Perturbative QCD methods have been previ- 
ously shown |ll[] to be a robust framework for describing 
these processes and earlier work jl0| had already shown 
that from the dominant decay mechanism, Fig. (1), the 
observed rate for B — > K*j jij was consistent with a 
pQCD approach. Using this experience, the question of 
the relative importance of subdominant processes to the 
decays B — > K*^ and B — > pr/ have been herein ad- 
dressed at the leading twist, leading a s level. By work- 
ing within a single, coherent framework the possibility 
of double counting present when using a mixed parton- 
meson approach has been avoided. 
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In the case of B — > K*j, subdominant charm quark 
loops 0, Mc E i 00 p, here called double-penguins (and else- 
where in the literature as "long-distance", vector 
meson like states) have been found to yield an enhance- 
ment of up to 30% in the total decay rate, depend- 
ing upon bound-state parameters of the K*. The un- 
certainty is due primarily to phase differences between 
-^ypeng, the dominant mechanism, and M cS \ oop . 

In the case of B — > pj, which depends upon the un- 
known CKM matrix elements Vb u and Vtd, additional 
subdominant processes involving up quark loops and also 
annihilation diagrams |L9| may no longer be CKM sup- 
pressed. Assuming SUf(3) symmetry of the p and K* 
distribution amplitudes, one can estimate the accuracy 
to which one could extract || |Vt<j/V ts | using the ratio of 
branching rates r B +_ >p + y /T B +^ K ,+ _,. A conservative 
estimate using present bounds Q for Vb c V* d /VbtV t * d and 
reflecting the sensitivity of this analysis to the various pa- 
rameters entering the calculation suggests that one can 
extract the ratio |Vtd/Vt s | to within a factor of two. 
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FIG. 1. The leading contribution to the decay B — > if* 7. 
The wavy line is a photon, the curly line is a gluon. 



FIG. 4. The function f(r = —2.25) (bottom curves) and 
f(r = —.5) (top curves) for the ratio F B +^ p + 7 /r s +^ K *+ 7 
as a function of Aqcd ■ The full lines are the results using the 
asymptotic distribution amplitude for the vector mesons; the 
dotted lines are those using the CZZ distribution. 



FIG. 2. One of the two "double-penguin" contributions to 
the decay B — > if* 7. The second (or photon gluon crossed) 
graph is required for gauge-invariance and to obtain a ultra- 
violet finite result. 



FIG. 3. The leading annihilation graph. The photon is 
emitted from the light quark in the B-meson. 



FIG. 5. f(r = -2.25) and f(r = -.5) for 
T B +^ 7 /r s +^ Jf ,o 7 . Notation same as in Fig. (4). 
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